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Inelastic neutron scattering on deuterated single-crystal samples is used to study Haldane-gap 
excitations in the new 5 = 1 one-dimensional quantum antiferromagnet Ni(C5Hi4N2)2N3(PF6) 
(NDMAP), that was recently recognized as an ideal model system for high-field studies. The Haldane 
gap energies A^, = 0.42 meV, Ay = 0.52 meV and A^ = 1.86 meV, for excitations polarized along the 
a, b, and c crystallographic axes, respectively, are directly measured. The dispersion perpendicular to 
the chain axis c is studied, and extremely weak inter-chain coupling constants Jy — 1.8 ■ 10~'' meV 
and Jx = 3.5 • 10~* meV, along the a and 6 axes, respectively, are determined. The results are 
discussed in the context of future experiments in high magnetic fields. 
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I. INTRODUCTION 

The unique properties of the one-dimensional (ID) 
integer-spin Heisenberg antiferromagnet (HAF) have 
captivated the minds of condensed matter physicists for 
the last 20 years. In total defiance of the quasi-classical 
picture of magnetism, the ground state in this system is 
a spin-singlet- a "quantum spin liquid" with only short- 
range (exponentially decaying) spin correlations. The 
excitation spectrum is also rather unique, and, even in 
the absense of any magnetic anisotropy, features a so- 
called Haldane energy gap.EJ The wealth of theoretical 
and experimental results accumulated to date provide a 
fairly complete understanding of the physics involved, 
and few mysteries remain, as far as the behavior of 
the idealized model is concerned. Particularly reveal- 
ing was the neutron scattering work done on real quasi- 
ID S" = 1 HAF compounds such as CsNiClaJa NENP 
(Ni(C2H8N2)2N02(C104)),B and Y2BaNi05 (Ref. |). 

One of the few remaining unresolved issues is the be- 
havior of Haldane-gap antiferromagnets in high magnetic 
fields. An externaUBagnetic field splits the triplet of Hal- 
dane excitations,Q'0i3"la driving one of the modes to zero 
energy at some critical field He. At this soft-mode quan- 
tum phase transition, static long-range antiferromagnetic 
correlations appear. In essence, the external field sup- 
presses zero-point quantum spin fluctuations that are re- 
sponsible for the destruction of long-range order in the 
spin-liquid state, and restores long-range spin coherency. 
Despite the success of high-field bulk studies of several 
material, it has been highly frustrating, that for most 
existing systems the actual values of the critical fields 
are unobtainable in neutron scattering experiments. As 
a result, the properties of the high- field phase are still 
largely unknown. It is even unclear, whether or not the 
high-field phase is commensurate, and virtually nothing 
is known about the excitation spectrum. 




FIG. 1. A schematic view of the antiferromagnetic spin 
chains in the NDMAP crystal structure. Only Ni, N and 
C atoms are shown. 

Traditionally, NENP was the workhorse of high-field 
studies. For this compound large single crystal samples 
can be fabricated, and the critical field He — 9 T (when 
applied parallel to the chain-axis) is, in principle, acces- 
sible in a neutron-scattering experiment. In practice, it 
is however easier to perform measurements with a mag- 
netic field applied perpendicular to the spin chains. The 
corresponding critical field for NENP is, unfortunately, 
much larger: He ~ 11 — 13 T. In any case, measure- 
ments deep within the high-field phase are not possible. 
Moreover, due to certain structural features, the phase 
transition in NENP is smeared out. The g-tensors of the 
S — 1 Ni^+ ions are staggered in this compound and lead 
to an effective staggered field when an external uniform 
field is applied.t2rt3 As a result, a static staggered mag- 
netization appears at arbitrary weak applied fields, and. 
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strictly speaking, no additional spontaneous symmetry 
breaking occurs at He- 

The recent j— . discov- 

ery of NDMAP (Ni(C5Hi4N2)2N3(PF6)),y a relatively 
easily crystallized Haldaiie-gap compound with a critical 
field of only around 4 T,t3 and no staggering of g-tensors 
within spin chains, promises to make the high-field phase 
readily accessible in neutron scattering studies. The crys- 
tal structure of this material is similar to that of NENP, 
and is visualized in Fig. |^. The AF spin chains run 
along the c axis of the tetragonal structure (space group 
Pnmn, a = 18.046 A, & = 8.705 A, and c = 6.139 A), 
and are composed of octahedrally-coordinated 5=1 
Ni^+ ions bridged by triplets of nitrogen atoms. These 
long nitrogen bridges account for a realtively small in- 
chain AF exchange constant J= 2.6 meV, as estimated 
from bulk x{T) measurements.Ej ESR and specific heat 
studies revealed a transition to the high-field phase at 
i/| = 3.4 T, and = 5.8 T, extrapolated to T ^ 0, 
for a magnetic field applied parallel and perpendicular to 
the chain-axis, respectivelyJlil'tj This anisotropy of crit- 
ical field is attributed to single-ion easy-plane magnetic 
anisotropy of type DS^. The anisotropy constant was 
obtained from bulk susceptibility data: D/J w 0.3. Us- 
ing the well-known numerical reSjult « 0.41J -|- 2pD^ 
« 0.41J - pD, p w 2/3j&B one can thus esti- 
mate the Haldane gap energies: A|| w 2.1 meV and 
A_L « 0.54 meV. 

While NDMAP appears to be an ideal model sys- 
tem for neutron scattering experiments in the high- 
field phase, additional characterization, particularly a 
measurement of inter-chain interactions and in-plane 
anisotropy, is required before such a study can be car- 
ried out. Finding the 3D AF zone-center is especially 
important, since it is at this wave vector that static long- 
range correlations are expected to appear in the high- 
field phase. Obviously, inelastic neutron scattering is the 
most direct method to obtain this information. In the 
present paper we report the results of a zero-field inelas- 
tic cold-neutron scattering study of deuterated NDMAP 
single-crystal samples, aimed at extracting this informa- 
tion. 



was roughly 0.3° and increased to 1.5° and 3° after the 
first and second cooling, respectively. Inelastic neutron 
scattering measurements were performed on the SPINS 
3-axis spectrometer installed on the cold neutron source 
at the National Institute of Standards and Technology 
Center for Neutron Research. Pyrolitic graphite crystals 
set for their (002) reflection were used as monochroma- 
tor and analyzer. Beam divergences were approximately 
40' - 80' - 80' - 240' through the instrument, with a 
cooled Be filter between the monochromator and sample. 
The measurements were done with a fixed final neutron 
energy Ef = 2.8 meV. The crystal was mounted with 
either the a or & crystallographic axis vertical, making 
(0, k, I) and {h, 0, 1) reciprocal-space planes accessible to 
measurements. The sample was cooled to T = 1.4 K in 
an " ILL-Orange" cryostat with a 70 mm diameter sample 
well. 
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II. EXPERIMENTAL PROCEDURES AND 
RESULTS 

Fully deuterated NDMAP single crystal samples were 
grown from solution as described in Ref. |l^. It was ob- 
served that the crystals tend to shatter when cooled to 
low temperature, and even more so when warmed back 
up. When wrapped in aluminum foil the crystals do not 
fall apart, but the width of the mosaic distribution in- 
creases dramatically with thermal cycling. Most inelas- 
tic neutron scattering data were collected on a 140 mg 
single crystal sample that was taken trough the cooling 
cycle only twice. The mosaic of the as-grown sample 



FIG. 2. Typical constant-Q scans measured at T = 1.2 K in 
an NDMAP deuterated single-crystal. Solid lines are based on 
global fits to the data, as described in the text. The shaded ar- 
eas represent partial contributions of the a-axis and b-axis po- 
larized Haldane gap excitations to the intensity. The dashed 
line indicates the background level. The magnetic excitations 
are resolution-limited and peak shapes are entirely defined by 
resolution effects. 

The primary goal of the experiment was to determine 
the a- and 6-axis dispersion relations and the polariza- 
tions of the two lower-energy Haldane-gap excitations, 
relevant for the soft-mode transition at He- Most of the 
data were collected in constant-Q scans at the ID AF 
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zone-center I = 0.5 in the energy transfer range 0-1 meV. 
In these dispersion measurements Q-resolution is partic- 
ularly important, so a flat analyzer was used. Typical 
scans for different momentum transfers perpendicular to 
the c-axis are shown in Fig. ||. At all wave vectors a well- 
defined peak corresponding to the lower-energy Haldane 
excitation doublet is clearly seen around 0.5 meV energy 
transfer. 
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FIG. 3. A constant-Qii scan measured in NDMAP using a 
horizontally-focusing analyzer configuration. Top axis shows 
the momentum transfer perpendicular to the chain-axis. 
Lines are as in Fig. |^. 

To observe the higher-energy Haldane-gap excitation 
and obtain an accurate measurement of the anisotropy 
constant we performed a constant-Q scan in the range 
0-2.4 meV (Fig. H). To maximize intensity we used a 
horizontally-focusing analyzer pointing c* towards the 
analyzer to maintain wave vector resolution along the 
chain. The scattering angle was varied so the projec- 
tion of wave vector transfer along the chain was 0.5c* 
throughout the scan. In addition to the peak seen in 
the low-energy scans, a weaker feature is observed at 
tku w 2 meV that can be attributed to the c-axis- 
polarized Haldane gap mode. 



III. DATA ANALYSIS AND DISCUSSION 

The measured constant-Q scans were analyzed using 
a parameterized model cross section, numerically convo- 
luted with the Cooper-Nathans 3-axis spectrometer res- 
olution function. Near the ID AF zone-center the single 
mode approximation (SMA) for the dynamic structure 
factor of isolated. Haldane spin chain is known to work 
extremely welLEIIIla For each channel of spin polarization 
the dynamic structure factor S^°'"'> {q,u!) can be written 
as: 



1 - cos(qc) Zv 
X d[nuj - nuja.q), 

(huja^qf ^Al+v^ sin2(q-c). 



(1) 
(2) 



Here v is the spin wave velocity, given hy v ~ 2.49J.li3 
The dimension-less constant Z defines the static stag- 
gered susceptibility of a Haldane spin chain: Xtt = 
Zw/A.^l In the above expression we have included the 
magnetic form factor f{q) for Ni^+ and the polariza- 
tion factor (^1 ~ ■ The structure factor for weakly- 
coupled chains can be |Calculated in the Random Phase 
Approximation (RPA)E3 The expression for 5("")(g>^) 
docs not change explicitly, but the excitations acquire 
dispersion perpendicular to the chain axis: 



{huJa 



■sm^{qc) + ZvJ'iq). 



(3) 



In this formula J'(q) is the Fourier transform of inter- 
chain exchange interactions, that we assume to be 
isotropic (of Heisenbetg type). According to numerical 
calculations Z w 1.26.EjThe form of J' can be guessed by 
looking at the crystal structure. The smallest inter-chain 
Ni-Ni distance (8. 705 A) is along the b crystallographic 
axis. We shall denote the corresponding exchange con- 
stant as Jy. The next-smallest inter-chain Ni-Ni distance 
(10.478A) is along the (0.5,0.5,0.5) direction. This in- 
teraction, however, is frustrated by in-chain AF interac- 
tions (any site in one chain couples to two consecutive 
sites in another chain), and is thus irrelevant, within the 
RPA, at momentum transfers qc « tt. Finally, the third- 
smallest inter-chain distance (18.046A) is along a. The 
corresponding coupling constant Jx is expected to be very 
small, due to the large "bond" length. The Fourier trans- 
form of intcr-chain interactions can thus be written as: 



J'{q) = 2Jx cos(qa) -I- 2Jy cos{qb). 



(4) 



The chain-axis exchange constant J is not very im- 
portant for our measurements: for scans collected at the 
ID AF zone-center it only slightly influences the peak 
shapes, due to a non-zero Q-resolution along the chain- 
axis. We can therefore safely use the value previously 
obtained from susceptibility measurements. The relevant 
parameters of the model are thus the three gap energies 
Aj,, Ay and A^, two inter-chain exchange constants Jx 
and Jy and an overall intensity prefactor. Three addi- 
tional parameters were used to describe the background: 
an energy-independent component, and the intensity and 
width for a Lorenzian profile centered at hto = to ac- 
count for incompletely resolved elastic scattering. 



*Here we have adopted the notation used in Ref. ^ 
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FIG. 4. Transverse dispersion of lower-energy Hal- 
dane-gap excitations in NDMAP, as determined in global fits 
to constant-<3 scans (solid lines), as described in the text. 
The data points represent excitation energies obtained in fits 
to individual scans. 

The difference between Aa; and Aj, , usually a result of 
in-plane magnetic anisotropy of type E{Sl — Sy), is too 
small to observe two and separate corresponding peaks 
in any single constant-Q scan. Distinct polarization fac- 
tors for the two branches however allow us to extract 
both parameters in a global fit to the data measured 
at different momentum transfers perpendicular to the 
chain-axis. As a first step, we simultaneously analyzed 
the data collected in (0, fc, Z) reciprocal-space plane (11 
scans with a total of 382 data points). A very good 
global fit using Eqs. |l| and ^ is obtained with a resid- 
ual — 1-6. The solid lines in Fig. |2^-c were calculated 
from the globally optimized parameter set. The shaded 
areas represent the contribution of each mode, and the 
dashed lines show the background level. Figure Q (right) 
shows the obtained 6-axis dispersion relation (solid line), 
that has a minimum at = 0.5. Symbols in this figure 
represent the excitation energies obtained in fits to in- 
dividual scans. A similar global fit (4 scans, 149 data 
points, = 1-7) was applied to all constant-Q scans 
measured in the (h, 0, 1) scattering plane (solid lines in 
Fig. |2|i and e). In this case the excitation energies at 
the zone-boundary (0, 0, 0.5) were fixed to the values ob- 
tained from the (0, fc, Z)-plane global fit. Only the a-axis 
exchange constant was refined. Dispersion of magnetic 
excitations along this direction is barely detectable. The 
dispersion relation obtained from our fits is nontheless 
plotted in a solid line in Fig. ^ (left). The fitting analysis 
suggests a shallow minimum at /i = 0.5. From Eq. ^ we 
can thus guess that the global minimum of the 3D dis- 
persion (3D AF zone-center) is located at (0.5,0.5,0.5). 
The energy gap for c-axis polarized excitations was de- 
termined by fitting the model cross section to the wide- 



range constant-Q scan measured with the horizontally- 
focusing configuration (Fig. ^, solid line). In this proce- 
dure the values of A^ and Ay were fixed. The param- 
eters determined by the analysis described above are as 
follows: A^ = 0.42(0.03) meV, Aj, = 0.52(0.06) meV, 
A, ^ 1.86(0.1) meV, Jy = 1.8(0.4) • lO^^ ^^y^ ^nd 
= 3.5(3.0) • lO-'i meV. 

The in-chain exchange constant J and anisotropy pa- 
rameter pQ_|ire readily obtained from the measured gap 
energiesBEa J « 0.81(A^ + + ^z) = 2.28 meV, 
D « 1{2A^ - A^ - Ay) = 0.70 meV, and D/J = 0.30, 
which is consistent with the bulk susceptibility result of 
Ref. [13 The relative strength of inter-chain interactions 
are Jy/ J « 7 • 10~^ and Jx/ J ~ 1.3 • 10~p. These ratios 
are very similar to those found in NENP.B 



IV. SUMMARY 

The small value of critical fields in NDMAP, com- 
pared to those in NENP, are a result of smaller in- 
chain exchange interactions and somewhat larger easy- 
plane magnetic anisotropy. The lowest-energy excita- 
tion in NDMAP is polarized along the crystallograhic 
a-axis. The 3D magnetic zone-center appears to be 
(0.5,0.5,0.5). Future experimental studies of the high- 
field phase should thus concentrate on this region of re- 
ciprocal space. 
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